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ABSTRACT 

In this paper, we use high-resolution smoothed particle hydrodynamics (SPH) sim- 
ulations to investigate the response of a marginally stable self-gravitating protostellar 
disc to a close parabolic encounter with a companion discless star. Our main aim is to 
test whether close brown dwarfs or massive planets can form out of the fragmentation 
of such discs. We follow the thermal evolution of the disc by including the effects of 
heating due to compression and shocks and a simple prescription for cooling and find 
results that contrast with previous isothermal simulations. In the present case we find 
that fragmentation is inhibited by the interaction, due to the strong effect of tidal 
heating, which results in a strong stabilization of the disc. A similar behaviour was 
also previously observed in other simulations involving discs in binary systems. As in 
the case of isolated discs, it appears that the condition for fragmentation ultimately 
depends on the cooling rate. 

Key words: accretion, accretion discs - gravitation - instabilities - stars: formation 
- planets: formation - brown dwarfs: formation 



1 INTRODUCTION 

Young stellar clusters are dynamic environments and en- 
counters between cluster me mbers can be quite co mmon 
(see, for example, Figure 6 of IScallv fc Clarkdl200lft . Such 
encounters can have a significant effect on the structure 
and e volution of the g aseous discs that surround young 
stars iBate et alJ 120031 . Some of these effects include a 
tidal truncation or dis ruption of the disc ijClarke fc Pringjel 
119931: lHall et alJll996ft. a burst of accre tion activity onto 
the central star jBonnell fc Bastienlll992T) . or the triggering 
of a gravitational instabili ty in the disc, which might then 
lead to disc fragmentation feoffin et alll998l : IWatkins et all 
Il998allbl) . The latter possibility has been often invoked as 
a possible formation mechanism for low mass companions, 
such as bro wn dwarfs, or massive p l anets (for some recen t 
reviews, see IWhitworth et all 120061 : iGoodwin eJ L alj j2006jh 
However, previous analys es of this process IIBoffin et all 
ll99Sl:IWatkins et alJll998allbh were limited by low resolution 
(which did not permit the proper resolution of the disc verti- 
cal structure) and by the replacement of the energy equation 
by a simple barotropic equation of state. All these analyses, 
in fact, assumed that the disc is isothermal, which might be 
a reasonable assumption at very large distances (the discs 
considered in these works are ~ 1000 AU in size), but is 
certainly not adequate for smaller discs. Indeed, as clearly 



stated in lWhitworth et alJ i2006h "it is important that such 
simulations be repeated, with a proper treatment of the en- 
ergy equation... to check whether low-mass companions can 
form at closer periastra". This is precisely the aim of the 
present paper. 

The importance of cooling in determining the condi- 
tions under which a gravitationall y unstabl e disc f ragments 
has been clearly pointed out bv iGammiel i200lT) initially 
from local shearing-box simulations and l ater confirmed in 
global disc simul ations (iRice et alJliool iRice et alJliool : 
iMeiia et alJl2005|) . The general result is that fragmentation 
occurs only if the cooling timescale is smaller than a few 
times the dynamical timescale in the disc. This result can 
be understood in the following way. The linear stability of a 
rotating disc against axisymmetric gravitational perturba- 
tions is described by the well known Toomre parameter: 



where c s is the sound speed, « is the epicyclic frequency 
(equal to the angular velocity f2 in a Keplerian disc) and E 
is the disc surface density. The disc is unstable if Q < 1, 
so that local perturbations will grow on a timescale of the 
order of the dynamical timescale The linear stability 

then only depends on the disc temperature (through c a ) and 
density E. However, the non-linear growth of the perturba- 
tions heats up the equilibrium state (through compression 
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and shocks), increasing the value of Q and hence stabilizing 
the disc, and preventing fragmentation, unless the heat pro- 
duced by compression can be removed efficiently. Since the 
perturbation grows on the dynamical timescale, if we want 
to have fragmentation, we require that cooling acts on the 
same timescale. 

Actually, an alternative way to destabilize the disc is 
to increase its surface density (rather than decreasing its 
temperature), which might result from the tid al interaction 
with a co mpanion star. In a ser ies of papers, iBoffin et alJ 
l|l99d) and lWatkins et all Jl998afbh have explored this pos- 
sibility by simulating isothermal discs which are initially sta- 
ble (with Q > 1 throughout the disc) and are destabilzed 
by stellar encounters. The process appears to be successful 
in that in general the disc can be effectively destabilized. 
However, in this case, the subsequent evolution of the sys- 
tem is essentially dictated by the isothermal condition: the 
perturbation will grow and produce substellar fragments. 
This is indeed what is observed in these simulations. How- 
ever, it now becomes essential to determine whether - in 
the presence of heating and cooling - the strong, non-linear 
perturbation induced by the encounter would still be able 
to produce fragmentation, or whether, similarly to the case 
where the perturbation is induced more gently, through a 
slow cooling, fragmentation is instead inhibited. 

Similar issues have also been considered, often includ- 
ing the effects of tidal heating and radiative cooling, in the 
context of planet formation in binary systems. However, the 
results up to now appear to be c ontradictory. Indeed , if on 
the one hand iNelsonl J2000T) and iMaver et at! (1200511 claim 
that the presence of a co mpanion inhi bits the formation of 
planets via fragmentation, |BossJ (|2006j), on the other hand, 
argues in favour of fragmentation. 

In this paper, we present the results of a series of 
smoothed particle hydrodynamics (SPH) simulations of the 
encounter of a star surrounded by a gaseous disc with a dis- 
cless companion. The gas disc is here allowed to heat up 
through compression and shocks and to cool down slowly, 
so that in isolation the disc does not fragment but attains 
a marginally stable, self-regulated state. We explore param- 
eter space by varying the orbital properties of the stars, 
their mass ratio, and the mass ratio between the star and 
the disc. None of our simulations resulted in fragmentation, 
thus proving that also in this situation the main require- 
ment to form substellar companions is that the cooling be 
fast. The disc is strongly perturbed and its density can be 
enhanced significantly, but the heat provided by compres- 
sion strongly overcomes the density increase, resulting in a 
substantial stabilization of the disc. 

The paper is organized as follows. In Section |2] we de- 
scribe the SPH code that we use and the physical setup of 
our simulations. In Section[3]we describe the evolution of our 
reference run. In Section 0] we discuss the effects of changing 
the main parameters of the problem. In Section|S]we discuss 
our results and draw our conclusions. 



2 NUMERICAL SETUP 
2.1 The SPH code 

Our three-dimensional numerical simulations are carried out 
using SPH, a Lagrangian hydrodynamic scheme 



lMonagharll992l). The ge neral i mple mentation is very similar 
to lLodato fe Ricd <2004 12005ft and iRice et all <2005l) . The 

gas disc is modelled with 250,000 SPH particles (500,000 in a 
run used as a convergence test) and the local fluid properties 
are computed by averaging over the neighbouring particles. 
The central star and the perturber are modelled as point 
masses onto which gas particles can accrete if they get closer 
than the accretion radius, taken to be equal to 0.5 code units 
for the central star and the high mass perturber (see below) 
and 0.25 code units for the small mass perturber. 

The gas disc can heat up due to pdV work and artificial 
viscosity. The ratio of specific heats is 7 = 5/3. Cooling is 
here implemented in a simplified way, i.e. by parameterizing 
the cooling rate in terms of a cooling timescale: 



dui 
~dt 



Ui 

"7 ' 

'-cool 
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where u; is the internal energy of a particle and the cooling 
timescale t coo i is assumed to be proportional to the dynam- 
ical timescale, t coo i = 0Q -1 . W e know from previous work 
jGammidl200ll ; IRice et aljl2005h that if < 3 - 7 (depend- 
ing on the ratio of the specific heats) the disc is unstable to 
fragmentation. Our aim is to establish whether a disc that 
would not fragment in isolation can be driven to fragmen- 
tation as a result of a stellar encounter. We have therefore 
set = 7.5, so that in isolation our disc is not expected to 
fragment. 

2.2 Disc setup 

The main physical properties of the disc at t he beginning of 
the sim ulat ion are again similar to those of lLodato fc Ricel 
J2004I) and lLodato fc Ricd J2005I) . The disc surface density 
E is initially proportional to _R _1 (where R is the cylindri- 
cal radius), while the temperature is initially proportional 
to R~ 1/2 . Given our simplified form of the cooling function, 
the computations described here are essentially scale free 
and can be rescaled to different disc sizes and masses. For 
reference, we will assume that the unit mass (which is the 
mass of the central star) is 1M@ and that the unit radius is 
1AU. In these units the disc extends from Ri n = 0.25AU 
to Rout = 25AU. The normalization of the surface den- 
sity is generally chosen such as to have a total disc mass 
of Mdisc = O.IMq (in a couple of cases, however, we have 
also considered different disc masses, see below), while the 
temperature normalization is chosen so as to have a mini- 
mum value of Q = 2, which is attained at the outer edge of 
the disc. This is obtained by setting T ~ 45 K at the outer 
edge. 

Initially, the disc is evolved in isolation for several outer 
dynamical timescales. The gen eral features of this initial 
evolution is described in detail in lLodato fc Rice) i2004f) . The 
disc starts cooling down until the vertical scale-length H is 
reduced such that H/R » M disc /M* = 0.1. At this point the 
disc becomes Toomre unstable and develops a spiral struc- 
ture that heats up the disc and maintains it close to marginal 
stability, i.e. close to Q = 1, roughly independent of radius. 

2.3 Adding the perturber 

Once the disc has reached the quasi-steady, self-regulated 
state of marginal stability, we add the perturber star. This 
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Ma- IMrs 
iw disc/ 1V1 Q 


No §3-s pcLrticlcs 


il,J pcrt / ivl 


Orbit 


Direction 


eccentricity 


r -/ATI 

' peri/ rtu 


SI 


0.1 


250,000 


0.1 


coplanar 


prograde 


1 


17 


S2 


0.1 


250,000 


1.0 


coplanar 


prograde 


1 


17 


S3 


0.1 


250,000 


0.1 


coplanar 


retrograde 


1 


17 


S4 


0.1 


250,000 


0.1 


'head on' 




1 


17 


S5 


0.1 


250,000 


0.1 


coplanar 


prograde 


7 


17 


S6 


0.05 


250,000 


0.1 


coplanar 


prograde 


1 


17 


S7 


0.2 


250,000 


0.1 


coplanar 


prograde 


1 


17 


S8 


0.1 


250,000 


0.1 


coplanar 


prograde 


1 


30 


S9 


0.1 


250,000 


0.1 


coplanar 


prograde 


1 


40 


HI 


0.1 


500,000 


0.1 


coplanar 


prograde 


1 


17 



Table 1. Summary of the main physical properties of the simulations carried out. 



is initially put at a distance of from the primary 

star, and is placed in a parabolic orbit. At this distance, the 
sudden addition of the perturber does not cause any signifi- 
cant perturbation to the disc structure. We have performed 
a large number of simulations, by varying the orbital param- 
eters, considering coplanar encounters (both prograde and 
retrograde with respect to the disc rotation) and 'head on' 
encounters, where the orbital plane of the secondary is per- 
pendicular to the disc plane. We have also performed a sim- 
ulation of a hyperbolic (coplanar) encounter. We have varied 
the disc mass, the perturber mass M pcrt and the pericenter 
of the perturber orbit r por i and we have also run simulations 
with different numerical resolutions. TableQsummarizes the 
main physical setup of the various simulations that we have 
done. 



2.4 Resolution issues 

One important aspect that needs to be taken into account 
is whether our simulations have enough resolution to re- 
solve the fragments, if fragmentation occurs. As shown by 
iBate fc BurkertJ jl997f) . SPH correctly reproduces fragmen- 
tation if the Jeans mass is larger than the minimum resov- 
able mass, defined as: 



M min = 2JV nciKh mi = 2M U 



^ ticigh 
JVtot 



(3) 



where mi is the mass of a single SPH particle, M to t (= Mdi sc 
in this case) is the total mass of the gas being simulated, iVtot 
is the total number of SPH particles used, and N ne i g h is the 
number of SPH neighbours (which is always kept close to 
50). 

To estimate the expected Jeans mass Mj, we use the 
fact that the most unstable wavelength to gravitational in- 
stability is of the order of the disc thickness H . We then 
have Mj « EH 2 « M diBC (H/R) 2 , from which we get: 



Mj 
M min 



N u 



N n 



1 /Mdis 

2 V M* 



neigh 



(4) 



For most of our simulations, we use 250,000 particles, 
which would then lead to Mj « 25M m j n . We thus have more 
than enough resolution to properly address the issue of frag- 
mentation. However, to be sure, we have also performed a 
convergence test of our results, by running a simulation with 
twice as many particles. The results of this convergence test 
are described below in Section 2] 



The second aspect related to resolution is that we re- 
quire artificial viscosity to play a role only when modeling 
shocks. In order to ensure this, we then require that the 
velocity difference accross a smoothing kernel is subsonic, 
i.e. hSl < Cs, where h is the smoothing length. This in turn 
requires that the smoothing length is smaller than the disc 
thickness H = c s /il. We have indeed checked that, even at 
the lower resolution of 250,000 particles, the average smooth- 
ing length is a fraction as 0.4 of the disc thickness. 



3 THE REFERENCE RUN 

In this section we describe the results of our reference sim- 
ulation, called SI, for which the main parameters are listed 
in the first line in Table Q Fig. Q shows images of the disc 
structure at different times during the evolution, where t = 
indicates the time at which the perturber reaches pericen- 
ter. The colour scale indicates the logarithm of disc surface 
density between 1 g/cm 2 and 10 4 g/cm 2 . The linear size 
of the images is 160 AU and the images are centered on 
the instantaneous position of the primary. The two stars are 
visible as white dots in this image. The perturbing star is 
initially in the top right hand corner of panel a. It can be 
seen that the disc starts off in a marginally stable state and 
displays a clearly visible spiral structure at t ~ —90.7 years 
(panel a). The disc orbits in an anti-clockwise motion. At 
t ~ —22.3 (panel b), the perturber begins to interact with 
the disc. At t — (panel c), the perturber reaches the peri- 
centre. After pericenter passage, the disc structure is signif- 
icantly perturbed, as shown at t ~ 19.1 years (panel d). At 
t ~ 81.2 years (panel e), there is barely any evidence of the 
pre-encounter spiral structure and the only visible structure 
is a pronounced tidal tail. The perturber has also carried 
away roughly 4% of the primary disc mass, most of it be- 
ing accreted on the perturber and the remainder forming 
a small circum-secondary disc, as well as causing some gas 
(~ 6 % of the disc mass) to be scattered away from both 
stars. At t — 302.4 years (panel f), the spiral structure is 
completely gone leaving a highly stable larger disc in which 
fragmentation has been prevented. 

Figure [5] shows the radial profile of Q for the refer- 
ence case for three different times: the start (solid) and end 
(dashed) of the simulation, and the time when the perturber 
reaches the pericentre (dotted). It can be seen that when the 
perturber reaches the pericentre it has stabilized roughly the 
whole outer disc. By the end of the simulation, the disc has 
become stable at all radii. Figure [3] shows the corresponding 




Figure 1. Images of the reference simulation at (from top left to bottom right) (a) t ~ —90.7 yrs (b) t ~ —22.3 yrs (c) t = yrs (d) 
t ~ 19.1 yrs (e) t ~ 81.2 yrs (f) t ~ 302.4 yrs 
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Figure 2. Radial profile of Q at t ~ —90.7 yrs (solid line), t ~ 
yrs (dotted line) and r. = 302.4 yrs (dashed line). The vertical 
dashed line shows the pericentre distance. 
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Figure 3. Radial profile of E at t ~ —90.7 yrs (solid line), f. ~ 
yrs (dotted line) and t = 302.4 yrs (dashed line). The vertical 
dashed line shows the pericentre distance. 



plots for the surface density, E. There is not much evolu- 
tion in the surface density before the perturber reaches the 
pericenter. However, it can be seen that E has become sig- 
nificantly steeper at the end of the simulation, suggesting 
an increase in the accretion rate. 

A more detailed view of the effects of the interaction 
on the disc structure can be obtained from Fig. 2] (solid 
line), where we plot the evolution of the sound speed and 
Q at a radius R r — 8.5 AU, i.e. halfway between the peri- 
center and the primary star; the moment at which the per- 
turber reaches pericenter is marked with a dashed line. Sev- 
eral things can be noted from these plots. First of all, the 
interaction produces a strong increase in sound speed, with 
the local temperature increasing from about 10 K to more 
than 100 K within a few decades of pericentre. Thereafter 
the sound speed declines on the local thermal timescale. The 
right hand panel shows that Q tracks the sound speed al- 
most exactly, which is consistent with the fact that changes 
in E at this radius are minor (a few tens of per cent at 
most). Evidently, the interaction has increased the specific 
entropy of the material, which must have been mediated by 
shock heating. Clearly, the net effect of the interaction is to 
increase Q and thus to stabilise the disc. 

4 EXPLORING THE PARAMETER SPACE 

4.1 Changing the direction of the perturber's 
orbit 

Figure [I] shows a comparison of sound speed and Q at 
R — 8.5AU for the prograde, retrograde and non-coplanar 
cases (simulations, SI, S3, and S4 respectively). The non- 
coplanar orbit is set up in such a way so that it impacts 
the disc at 90°. It can be seen that the largest effect on the 



disc structure is obtained with a prograde passage, while 
the retrograde and especially the 'head on' encounter have 
a substantially reduced effect. The fact that the 'head on' en- 
counter has a smaller effect can be easily understood since 
the timescale over which the perturber interacts with the 
disc is significantly reduced (a similar effect occurs also in 
the case of a hyperbolic encounter, see below). The com- 
parison between the two coplanar encounters (pr ograde and 
retrog rade) is in agreement with the results of lHall et alJ 
(1996), who have shown that retrograde encounters induce 
a smaller amount of angular momentum transfer between 
the perturber and the disc. We also note that the timescale 
on which the temperature reaches its peak is roughly double 
that for the prograde encounter. 

4.2 Changing the disc mass 

Figure |S] compares the results of simulations SI, S6 and S7, 
in which we change the disc mass. The sound speed profiles 
show something interesting: even though the temperature of 
the higher mass disc is overall the largest, its increase with 
respect to the pre-encounter value is the lowest, while it is 
the largest for the lighter disc. This also shows up in the 
plots of Q, where for the high disc mass case, the increase 
in Q after the encounter is smaller than for the lower disc 
mas ses. 

iMaver et al.l (l2005t) . in their analysis of gravitational in- 
stability in binary systems suggested that a lighter disc is 
more likely to fragment when perturbed by a binary com- 
panion as the heating caused by shocks is not enough to bal- 
ance the large increase in density, thus resulting in a less sta- 
ble disc. However, the present simulations show that a lower 
disc mass is stabilized more efficiently than a heavier disc. 
It appears that the smaller fractional temperature increase 



6 Lodato, Meru, Clarke & Rice 




Figure 4. Azimuthal averages at _R = 8.5AU of sound speed (left) and log(Q) (right) as a function of time for the simulations involving 
prograde (solid line), retrograde (dotted line) and non-coplanar (dashed line) orbits. 
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Figure 5. Azimuthal averages at R = 8.5AU of sound speed (left) and log(Q) (right) as a function of time for simulations SI (solid line; 
M disc = O.lAf©), S7 (dashed line; M disc = 0.05M Q ) and S6 (dotted lne; M disc = 0.2M Q ). 



for a high disc mass, with comparable effect on E, makes 
a higher disc mass relatively more susceptible to fragmen- 
tation. Nevertheless, we emphasise that for all disc masses 
studied, Q increases during the encounter, i.e. the effect of 
the encounter is one of stabilisation. 



is indeed confirmed in Fig. [7] that shows the time evolution 
of the disc surface density for the reference case (solid line) 
and for simulation S2 (dotted line). While in the reference 
case after the encounter the surface density is roughly the 
same as the pre-encounter value, in the high perturber mass 
case, it is decreased by roughly a factor 3. 



4.3 Changing the perturber's mass 

Figure |S| compares the Reference simulation with the simu- 
lation with a perturbing star that is ten times heavier (sim- 
ulations Sf and S2). As expected, increasing the mass of the 
perturber has a much larger effect on Q, and sound speed, 
consistent with more vigorous shock heating. We note that 
Q does not in this case just track the evolution of the sound 
speed, implying that the effect of the higher mass is to de- 
crease the local column density by a factor of a few. This 



4.4 A hyperbolic encounter 

In Fig. |H]we show the results of a hyperbolic encounter (sim- 
ulation S5). In general, the effect is similar to the parabolic 
one, but smaller in amplitude. This is in general consistent 
with the fact that the interaction occurs on a much faster 
timescale for a hyperbolic encounter and is therefore less 
effective in perturbing the disc structure. 
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Figure 6. Azimuthal averages at R = 8.5AU of sound speed (left) and log(Q) (right) as a function of time for simulations SI (solid line; 
A/p Cr t = 0.1M©) and S2 (dotted line; M peIt = M©). 




Figure 7. Azimuthal averages at R = 8.5AU of surface density £ 
as a function of time for simulations SI (solid line) and S2 (dotted 
line; high perturber mass). 



4.5 Changing the pericenter distance 

Figure [3] shows that overall, encounters with larger pericen- 
ter distances have a much smaller effect on the temperature 
and Q. There is slight increase in temperature for the case 
where -R per i = 30, whereas for the 7? pG ri = 40 case also the 



temperature is barely modified. 



4.6 A convergence test 

Finally, in Fig. 1101 we show the results of our convergence 
test. We have re-run simulation SI at a higher resolution, 
by doubling the number of SPH particles used. Figure 1101 
clearly shows that there is little effect on the temperature 
and the Toomre parameter when the number of particles is 
increased to 500,000. Therefore, we can conclude that the 



resolution of the reference simulation is adequate and does 
not affect the outcome of the results. 



5 DISCUSSION AND CONCLUSIONS 

In this paper, we have analysed the effects of an encounter 
with a discless companion on the structure and evolution of 
a protostellar disc. We have improved on previous analyses 
by considering in more detail the energy balance for the gas, 
thus going beyond the approximation of isothermal evolu- 
tion used in the past. This turns out to be very important, 
because the outcome of the encounter is strongly dependent 
on the tidal heating induced in the disc. Indeed, probably the 
main result of the present paper is the demonstration that an 
encounter with a companion inhibits rather than promotes 
fragmentation of a gravitationally unst able disc. This is i n 
marked contrast with previous results jBofEn et al.lll99ci) . 
who instead have shown effective fragmentation following 
an encounter . The main d i fferen ce between ou r simulations 
and t hose of iBoffin et alJ il99ST) is that while iBoifin et alJ 
( 1998) consider much larger discs, for which the isothermal 
approximation is probably adequate, we instead follow the 
heating and cooling processes in the disc and therefore allow 
the growth of gravitational instabilities to feed back upon 
the thermodynamic state of the disc. We include heating 
from pdV work and shocks and cool the disc down with a 
cooling rate that is sufficiently small that the disc would not 
fragment in isolation. 

We thus conclude that the main parameter that deter- 
mines the fragmentation of discs is the cooling rate, regard- 
less of the way in which the disc is driven to instability (ei- 
ther through cooling, as in the simu lations of lGammieH200ll : 
iMaver et al l2002HRice et alJl200l or by a dynamical inter- 
action, as simulated here). A disc that does not fragment in 
isolation, cannot be driven to fragmentation through an in- 
teraction with a companion. Indeed, if anything, tidal heat- 
ing could prevent fragmentation even in cases where the disc 
would otherwise fragment, as shown bv lMaveT et all d2004l) 
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Figure 8. Azimuthal averages at R = 8.5AU of sound speed (left) and log(Q) (right) as a function of time for simulations SI (solid line; 
parabolic encounter) and S5 (dotted line; hyperbolic encounter. 



(and previously, but in two dimensions, byjNelsonl^OOO ) , 
in the context of discs in a binary system. iBosall^'oOrj) . on 
the other hand, has shown that in binary systems, disc frag- 
mentation is indeed possible if the cooling timescale is suffi- 
cie ntly small ( estimated to be about 1 — 2 orbital timescales 
by|Bc«s!l200a). We note however, that the discs modelled 
bv lBosa pOOfi) were also shown to fragment even in the ab- 
sence of binary companions. Here again, it would therefore 
appear, that the cooling timescale is the main determinant 
of whether discs fragment or not and that the role of com- 
panions in either promoting or suppressing fragmentation is 
rather minor. 

Our simulations adopt a simplified prescription for the 
cooling. In reality, the cooling timescale is going to depend 
on disc properties, such as density and opacity. In order 
to determine the stability of discs against fragmentation it 
is therefore needed to carefully evaluate the actual cooling 
rate. In the context considered here , it is then interesting 
to note that IWhitworth et all i2006l) have considered this 
issue analytically and have shown that, with realistic opacity 
laws, the effect of an encounter is to reduce the cooling rate, 
further reinforcing the conclusions obtained in our work. 

The simulations presented here involve an encounter 
between a protostellar disc and a disc-less perturber. 
Our results are th erefore directly comparable to those of 
Boffi n et afll|l99tf) . It is however also interesting to consider 
the more complex behaviour and the much larger parame- 
ter space involved in d isc-disc interactions, as considered by 
IWatkins et all l)l998allb|) . We postpone this investigation to 
future analyses. 
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Figure 10. Azimuthal averages at R = 8.5AU of sound spped (left) and log(Q) (right) as a function of time for simulations SI (solid 
line; 250,000 particles) and HI (dotted line; 500,000 particles). 



